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We study the ordering of the spin and the chirality in the fully frustrated XY model on a square 
lattice by extensive Monte Carlo simulations. Our results indicate unambiguously that the spin and 
the chirality exhibit separate phase transitions at two distinct temperatures, i.e., the occurrence 
of the spin-chirality decoupling. The chirality exhibits a long-range order at T c = 0.45324(1) via 
a second-order phase transition, where the spin remains disordered with a finite correlation length 
£s(Tc) ~ 120. The critical properties of the chiral transition determined from a finite-size scaling 
analysis for large enough systems of linear size L > f s (T c ) are well compatible with the Ising 
universality. On the other hand, the spin exhibits a phase transition at a lower temperature T s = 
0.4418(5) into the quasi-long-range-ordered phase. We found T](T S ) — 0.201(1), suggesting that the 
universality of the spin transition is different from that of the conventional Kosterlitz-Thouless (KT) 
transition. 
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I. INTRODUCTION 

Fully frustrated XY (FFXY) model is a basic statis- 
tical mechanical model of geometrically frustrated sys- 
tems P, An example is the antiferromagnetic XY 
spin system on the two-dimensional (2D) triangular lat- 
tice [3j. Another realization is a Joscphson junction array 
(JJA) on the square lattice under an external magnetic 
field such that the flux density per plaquette is / = 1/2 
[3, [1[ . These systems possess an Ising- like Zi symmetry 
with respect to global spin reflections in spin space and 
a continuous U(l) symmetry with respect to global spin 
rotations in spin space. 

The ground state of the model is characterized by the 
opposite sense of the chirality. The chirality is an Ising- 
like variable taking either plus or minus corresponding 
to the right- and left-handed noncollinear spin structure 
[l[ . More precisely, in the ground state of the model, the 
chirality exhibits a chcckcr-board-likc pattern much as in 
the antiferromagnetic Ising model on the square lattice. 
Since the chirality changes its sign with respect to global 
spin reflections in spin space, the Z2 spin-reflection sym- 
metry is broken in the ground state of the model. Then, 
by an analogy to the 2D Ising model, it would be natural 
to expect that the ordering of the chirality takes place at 
a finite critical temperature T c . Concerning the contin- 
uous U(l) spin symmetry, although the Mcrmin- Wagner 
theorem inhibits the onset of a true spin long-range or- 
der at any finite temperature [f|, a finite-temperature 
transition into a quasi-long-range-ordered state like the 
Kosterlitz-Thouless (KT) transition Q is still expected 
to occur. 

A central issue in the studies of the FFXY model, since 
the seminal work by Teitel and Jayaprakash , has been 
to clarify how these two distinct types of orderings take 
place. An interesting possibility among others, which 
has attracted researchers for decades, is that even at the 
temperature where the chirality establishes a long-range 



order the spin (superconducting phase in case of the 
JJA) may remain disordered due to thermally excited, 
unbound vortices. Then, the orderings of the chirality 
and of the spin take place at two separate temperatures 
T = T c and T = T s such that T c > T s . We emphasize 
that the possibility of such a "spin-chirality decoupling" 
is a basic general problem in the studies of frustrated 
magnets including spin glasses Q. The universality of 
the two transitions is also an interesting issue. However, 
the separation between T c and T s in the present system is 
extremely small, if any, which makes the solution of this 
problem technically difficult in spite of intensive studies 
done for more than two decades [9l-[35|. 

Under such circumstances, the purpose of the present 
paper is to clarify the nature of the phase transitions 
of the 2D FFXY model by performing extensive Monte 
Carlo (MC) simulations for systems as large as up to 
L ~ O(10 3 ). As was recently suggested by Hasenbusch, 
Pelissetto and Vicari (HPV) [2lTl22T|. such large system 
sizes are needed to extract asymptotic crucial proper- 
ties of the system. Although our system size is com- 
parable to Hsenbusch et al, wc study various physi- 
cal quantities including those which were not studied in 
Refs.(2~il [22l| . Namely, in addition to the Binder param- 
eters, the correlation lengths and the helicity modulus 
studied in Refs. j2lT [22| . we also investigate the specific 
heat, the chiral susceptibility and the vorticity modu- 
lus [H, in order to shed further light on the order- 
ing of the model. Furthermore, Refs.fU [22[ restricted 
their study just to the transition point, whereas we study 
an explicit temperature dependence of various physical 
quantities in a finite temperature range around the tran- 
sition point (s). 

Then, we have found convincing evidence that the two 
phase transitions take place at mutually close but two 
distinct temperatures. We have also found that the crit- 
ical properties of the chiral transition are well compat- 
ible with the Ising universality class. We corroborate 
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Refs.(2ll. 1221 ] in these conclusions. On the other hand, 
we have found that the universality of the spin transi- 
tion is not compatible with the standard KT transition, 
a conclusion different from that of Refs.j2ll.l22l]. 

The organization of the paper is as follows. In the 
next section, we introduce our model and explain the 
method of our simulations together with various physical 
observables we calculate. In sec. IIII1 we present and 
discuss the results of our MC simulation. Finally in sec. 
IIV1 we conclude with a summary and some additional 
discussions. In appendix [A] some details of the coulomb- 
gas representation of the vorticity modulus are given. 



II. MODEL AND METHOD 

We study the fully frustrated XY (FFXY) model on 
a square lattice of size N = L x L. We label the sites of 
the square lattice as i = 1,2,. . . ,N and denote their x 
and y-coordinates as x% = 1, 2, . . . , L and j/j = 1, 2, . . . , L. 
Similarly, we label the plaqucttes as n = 1,2, ... ,N and 
denote their coordinates as x n = 1,2, ... ,L and y n = 
1,2,. ..,L. 

The Hamiltonian of the model is given by, 



We study the onset of the chiral order via the chiral 
Binder parameter, 

* = 5i X -5wJ' (4) 

and the chiral correlation length, 

p = 1 / (K(0)| 2 ) ' 

" c 2sin(7r/L)V (|m c (g min )| 2 ) ' (> 

where 

m M = jj52(-l) M e^K n . (6) 
n=l 

is the Fourier transform of the staggcred-chirality where 

Qmin ~ r L~&x &ud T n (^n?2/n)- 

In addition, we examine the staggered chiral suscepti- 
bility defined by, 

Xc = N/3(m 2 c }, (7) 
where j3 is the inverse temperature. 



H — ^2 JijSi ■ Sj 



(1) 



where Si = (Sf,S?) = (cos U sin 6$ (0 < 9 { < 2ir) is 
the two-component spin variable at the site i. The sum- 
mation (ij) is taken over all nearest-neighboring sites. 
The interaction bond takes two possible values = ± J 
such that every plaquette on the square lattice is frus- 
trated , i.e., the product of Jij around each plaquette 
is sign(Jjj) — — 1. In the following, we measure the 
(reduced) temperature k^T / J in a unit with k^/J = 1. 
We denote the thermal averages as (...). 



A. Chirality-related observables 

Let us introduce first physical observables associated 
with the chirality. We define the local chirality at a pla- 
quette n by, 



^= sign(Jij) sin(0i -0 ), 



(2) 



where the sum is taken over the four bonds around the 
plaquette. The normalization factor is chosen such that 
the absolute value of the chirality becomes unity in the 
ground state. The checker-board-like order of the chiral- 
ity can be detected though a 'staggcred-chirality', 



B. Spin-related observables 

Next, we introduce the physical quantities associated 
with the spin. The spin order parameter is defined via 
the sublattice magnetization. The square lattice is de- 
composed into four sublattices, and the sublattice mag- 
netization is defined by, 
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where a = 1, 2, 3, 4 denotes each sublattice. 

In 2D, the spin order parameter should vanish in the 
thermodynamic limit TV — > oo at any finite temperature. 
However, the spin may still establish a gwasi-long-range 
order below a finite transition temperature T s . To detect 
such a spin order, we measure the spin Binder parameter, 
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and the spin correlation length 

i r~< 



where 
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and ff m in = (2tt/L, 0). We also calculate the sublatticc 
spin susceptibility denned by, 



Xs 



(12) 



In addition, we study the helicity modulus and the 
vorticity modulus which detect the rigidity of the system 
against appropriate spin deformations. The two quanti- 
ties arc related but have physically different meanings as 
explained below. 

The helicity modulus T measures the rigidity of the 
system against an infinitesimal uniform twist of the spins 
A9 applied on the boundary. More precisely, it is defined 
through the change of the total free energy AF with re- 
spect to an infinitesimal twist on the spin configuration 
along, say, y-axis 0, — By — > 9{ — By + (A6/L) where i' is 
the nearest neighbour site of site i along the +y direction. 
One readily finds, 



AF = Y(A#) 2 + <9(A 4 
The helicity modulus is then given by, 
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On the other hand, the vorticity modulus [36|,[37|J mea- 
sures the rigidity of the system against a creation of an 
extra vortex. Disregarding the discrete nature of the vor- 
ticity, we consider an infinitesimal increase of the vor- 
ticity represented by the phase differences B{ — 6j — > 
9i — 9j + m<f>ij across each bond < ij > where m is as- 
sumed to be small. Here 4>ij = — 4>ji is the solid angle 
between site i and j with respect to an arbitrarily cho- 
sen origin in the system where the core of the vortex is 
supposed to be located. In practice, we define the solid 
angle as fyj — 



4>ij/v a 



Vi 



Here is the unit 



vector connecting the neighboring sites i and j, and 0y 
is a unit vector tangential to the circular path passing 
through the site i enclosing the origin. One readily finds 
that the free-energy increase can be expressed as, 



with 



V = 



AF = Vm 2 + 0(m 4 ), 



(15) 
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For large enough systems, we expect the scaling be- 
haviour, 



where v(T) is the vorticity modulus. In appendix \K\ 
we present a coulomb-gas representation of the vortic- 
ity modulus which might provide some insight into the 
physical meaning of this quantity. 

In practice, we evaluate the vorticity modulus v(T) in 
two different ways. The first way is to define the effective 
vorticity modulus by, 



vi(T,L) 



V(T,L) 
InL ' 



(18) 



and examine its L — >• oo limit. The second way is to 
define the effective vorticity modulus using the data at 
two different system sizes L\ and Li [37j . 



u 2 (T,ii,i 2 ) 



ln(Li/L 2 ) 



(19) 



V(T, L) w c(T) + v(T) InL, 



(17) 



and examine its L\, Li — > oo limit. The advantage of w 2 
compared with the first one v± (T) is that the anticipated 
constant term c(T) in Eq. P7|) is cancelled so that a 
finite-size correction becomes smaller. 



C. Monte Cairo simulations 

We use a Monte Cairo method which combines the 
Metropolis method and the over- relation method j38| . 
Within a single Monte Carlo Step (MCS), the spins 
are updated by one sweep using the Metropolis method 
which is followed by three consecutive over-relaxation 
sweeps. System sizes of L = 32, 64, 128, 256, 512, 1024 
plus L = 384, 768 under periodic boundary conditions 
are studied. We use the histogram method to evaluate 
physical quantities around the critical temperature. Typ- 
ical MC steps used for equilibration and measurements 
arc 10 8 MCS for L = 512, 768, 1024 around the chiral 
transition temperature and 5 x 10 6 MCS around the spin 
transition temperature. 



III. RESULTS 

In this section, we present the results of our MC sim- 
ulations. 



A. Specific heat 

We begin with the specific heat C = (3' 2 ((H 2 ) - (H) 2 ). 
As shown in Fig. QJa) > the specific heat exhibits a sharp 
peak as a function of the temperature at T = Tp ak(L). 
The height of the peak increases with increasing L, sug- 
gesting the occurrence of a second-order phase transition. 
Meanwhile, there is no indication of other phase transi- 
tions. This could be explained if one recalls that the 
expected spin transition might be of the KT-type which 
is known not to accompany a distinct anomaly in the 
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FIG. 1: The temperature dependence of the specific heat per 
spin. The panels a) and b) show the data around the chiral 
transition temperature T c and spin transition temperature T s , 
respectively. The arrows indicate the estimated transition 
temperature for the chirality T c , and for the spin T s . 



specific heat. Therefore it is natural to expect that the 
specific-heat-peak is associated with the chiral transition. 

In order to estimate the bulk chiral transition tem- 
perature T c , we show in Fig. [3] the size dependence of 
the specific-heat-peak temperature T pea k(X). Interest- 
ingly, the system-size dependence of the peak temper- 
ature Tpeak(i) exhibits a crossover behaviour around 
L = L x ~ 200: For sizes smaller than L x ~ 200, 
Tpeak(L) exhibits a near- linear dependence, while for 
sizes greater than L Xl T pea j £ (L) exhibits another near- 
linear dependence but with a different slope. By making 
a linear extrapolation of T pea k(£) for large enough sys- 
tem sizes L > L x ~ 200, we get an estimate of the chiral 
critical temperature, T c = 0.45324(1): See Fig. [3(b). It 
implies that the long-range order of the chirality is es- 
tablished at T c = 0.45324(1) via a second-order phase 
transition. Since the finite-size scaling ansatz implies 
Tp ea k(L) w T c + constL -1 /^ where the critical expo- 
nent v c describes the divergence of the chiral correlation 
length, the chiral transition is characterized by v c ~ 1, 
which is consistent with the the Ising universality class. 
Our present estimate of T c and the Ising nature of the 
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FIG. 2: Spin correlation length £ S (T,L) measured at the 
chiral transition temperature T c = 0.45324 for various system 
sizes. 



transition are consistent with the results of many recent 
works, particularly with the most recent result by Hascn- 
buseh et al who gave T c = 0.45324(1) [H|,[22j]. 

Now we turn to the issue of the crossover observed 
around L x ~ 200, which is likely to reflect a change in 
the nature of the ordering of the spin. In the inset of 
Fig. [2J we show the size dependence of the spin correla- 
tion length £ s at the chiral transition point T c — 0.45324 
determined above. As can be seen from the figure, £ s (7c) 
tends to saturate for large enough sizes of L > 200, 
yielding £ S (T C ) ~ 120, which is indeed comparable to 
the crossover length scale L x ~ 200 determined above. 
This result indicates that the observed size-crossover re- 
flects the change in the ordering behavior of the system 
from the 'spin-chirality coupling' regime where the spin 
correlation dominates the ordering to the 'spin-chirality 
decoupling' regime where the chiral correlation outgrows 
the spin correlation and dominates the ordering of the 
system. 

The universality of the chiral transition of the 2D 
FFXY model has been a matter of intense debate in 
the previous studies. Most probably, the source of the 
confusion is this size-crossover effect associated with the 
spin-chirality coupling/decoupling. If one looked at small 
enough systems compared with the crossover length L x , 
one should find an effective i/-exponent close to unity 
as demonstrated in Fig. [TJa), which would suggest the 
Ising universality. This was actually observed in simu- 
lations in early days [H, [24|, where system sizes in the 
range 10 < L < 50 were studied. However, if one focused 
on intermediate length scales around L x , one would find 
v < T which would suggest the non-Ising universality 
[H Qj, [13. Then, at length scales sufficiently 

larger than L x , the effective i/-exponent would become 
closer to unity again [U [22| . Our present estimates of T c 
and v are in perfect agreement with the recent estimates 
of Hasenbusch et al [2l|, [22| . 
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FIG. 3: The specific-heat-peak temperature T pea k(£) plotted 
versus the inverse linear size 1/L, for (a) all sizes investigated, 
and for (b) larger sizes. The straight lines represent linear 
fits of the data at smaller sizes L < 200, and at larger sizes 
L > 200. 



To get more insight into the ordering process of the 
model, we proceed to analyze in the next subsections 
physical quantities which are directly related to the order 
parameters. 



B. Chirality-related quantities 

In this subsection, we examine the ordering of the chi- 
rality. As shown in Fig. HJa), the chiral Binder parame- 
ter g c (T, L) of different sizes exhibit clear crossings. The 
scaled correlation-length, or the chiral correlation-length 
ratio £ C {T, L), also exhibits similar crossings as shown in 
Fig. Hb). 

To analyze the size dependence of the crossing points 
quantitatively, we label the system sizes in the ascend- 
ing order L\ < L 2 ■ ■ ■ and plot in Fig. [5] the crossings 
point T cross (L) of the Binder parameter and of the chiral 
correlation-length ratio at adjacent sizes and L n 

L n+1 +L„ 



as a function of the average size L av = " . As in 

the case of the spccific-hcat-pcak temperature, the cross- 
ing temperature T CIOSS (L) of both the Binder parameter 
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FIG. 4: The temperature dependence of (a) the chiral Binder 
parameter g c , and of (b) the chiral correlation-length ratio 
tc/L. 



and the correlation-length ratio exhibits a size-crossover 
at around L ~ 200, the same crossover length observed 
in the specific heat, where a non-monotonic size depen- 
dence with an abrupt turnover takes place. With use of 
the data of larger sizes L > L x , the chiral transition tem- 
perature is estimated as T c = 0.45324(5), which agrees 
well with the estimate above obtained from the specific 
heat T c = 0.45324(1). 

To examine the universality of the chiral transition 
further, we analyze the chiral susceptibility \c- At the 
chiral transition temperature T c , a power- law behaviour 
X C (T C , L) oc L 2 ~ Va is expected with a chiral anomalous- 
dimension exponent r/ c . We analyze the data at each 
temperature T shown in Fig. [6] by fitting them to the 
form Xc(T, L) = a(T) + b(T)L 2 - r >°«( T ') with three fitting 
parameters a(T), b(T) and i] c fi(T). 

The temperature dependence of the effective exponent 
rj e g(T) is displayed in Fig.[7][a). We find that, at the chi- 
ral transition temperature T c = 0.45324 estimated above, 
the effective exponent r] e g(T) takes a value 0.25(3). This 
observation is consistent with the expected Ising value 
r/ = 1/4, giving further support to the Ising nature of the 
chiral transition. 

In Fig. (7](b), we display the temperature dependence of 
our another fitting parameter a(T). Presumably, this pa- 
rameter reflects fluctuations at short length scales where 
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FIG. 5: Crossing temperatures of the chiral Binder parameter 
g c and of the chiral correlation-length ratio £,c/L are plotted 
versus the inverse linear size 1/L. The dashed line is a linear 
fit which yields T c = 0.45324(5). 
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FIG. 6: The chiral susceptibility divided by L 2 is plotted 
versus the system size L on a log-log plot. The solid line 
represents a fit to the form X c(T, L) = o(T) + &(T)L 2 -"<= ff (T) 
of the data at T = T c — 0.45324. For comparison, a power- 
law L~ v behavior with the Ising exponent r\ — 1/4 is also 
shown by the dotted line. 



chiral and spin orders are not separated. The increase of 
a(T) with decreasing temperature T may be interpreted 
as due to the growth of the spin correlation length £ B (T) 
approaching T s . 



C. Spin-related quantities 

Next, we examine the ordering of the spin. The 
spin Binder parameter g s (T, L) and the spin correlation- 
length ratio £ S (T,L)/L are shown in Fig. [8ja) and (b), 
respectively. Data of different sizes exhibit crossings as in 
the case of chirality shown in Fig. 0J From the raw data, 
it is already evident that the spin crossing temperatures 
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FIG. 7: The temperature dependence of the parameter val- 
ues obtained in the fitting analysis of the chiral susceptibility 
based on the form \o(T,L) = o(T) + b(T)L 2 -^ ii{ - T) : a) the 
effective exponent rj e ff(T), and b) the constant a(T). 



are considerably lower than chiral transition temperature 
T c = 0.45324(1). It means that the chirality and the spin 
exhibit two separate transitions. 

In addition, the raw data also suggests a qualitative 
difference in the nature of the orderings of the spin and 
the chirality. In case of the spin shown in Fig. |8l the sep- 
aration between the data points of different sizes tends 
to become increasingly small for larger systems, exhibit- 
ing a marginal behavior, whereas in case of the chirality 
shown in Fig. [4] the data points of different system sizes 
remain well separated even for larger systems. This dif- 
ference reflects the fact that, while the chirality exhibits 
a finite long-range order at finite temperature even in 2D, 
the spin cannot establish a true long-range order at finite 
temperature but only a quasi-long-range order. 

The effective vorticity modulus v\ and defined by 
Eq. (JTSJ) and Eq. ([T9j are displayed in Fig. g^a) and 
(b), respectively. The effective vorticity modulus v\ ex- 
hibits a stronger finite-size effect compared with V2, as 
expected. V2 exhibit merging behaviour at low tempera- 
tures suggesting the onset of the quasi-long ranged order. 
Again, it can be seen that the crossing temperatures of 
the vorticity modulus are considerably lower than the 
chiral transition temperature T c = 0.45324(1). 
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FIG. 8: The temperature dependence of (a) the spin Binder 
parameter g s , and of (b) the spin correlation- length ratio £ a /L. 
Evidently the crossing temperatures are lower than the chiral 
transition temperature T c = 0.45324. 



FIG. 9: The temperature dependence of the two types of the 
effective vorticity modulus, (a) vi defined by Eq. 1)180 . and b) 
V2 defined by Eq. (JTSj) . 



Now we estimate quantitatively the spin transition 
temperature T s based on our data of the crossing temper- 
atures of several observables presented above. Thus, the 
crossing temperatures T CIOSS (L) of the spin Binder param- 
eter g B (T, L), the spin correlation- length ratio £ a (T, L)/L 
and the effective vorticity modulus V\ between the two 
adjacent sizes L\ and Li are plotted against l/(lnL al ,) 2 
in FigJlOl Note that this form is motivated by the stan- 
dard KT form £ s ~ cxp[c/\/T — T s }. From the combined 
fit of these data where a common spin transition temper- 
ature T s is assumed, we get an estimate of the bulk spin 
transition temperature T s = 0.4418(5). 

In FigJTU] we also display the size dependence of 
the crossing points of the spin correlation-length ratio 
£ S (T, L)/L between the two sizes L and sL with s = 2, 4, 8 
plotted as a function of 1/ (In L) 2 . From the combined fit 
of these three data sets where a common spin transition 
temperature T s is assumed, we get T s = 0.4414(3) which 
is consistent with our estimate presented above. 

In order to examine the universality of the spin tran- 
sition, we also analyze the sublattice spin susceptibility 
Xs- At and below the spin transition temperature T s , 
a power-law behaviour Xs(T c , L) oc L 2 ~''( T ) is expected 
with a spin anomalous-dimension exponent rj(T). We 



analyze the data shown in Fig. [TT]by fitting them to the 
form Xs{T,L) = b{T)L 2 -^ tt{ - T ^ with two fitting parame- 
ters b(T) and rj e ff(T) at each temperature T. 

The temperature dependence of the effective expo- 
nent rj c s(T) is displayed in Fig. IT^T a). We find that, 
at the spin transition temperature T s = 0.4418(5) esti- 
mated above, the effective exponent r] c g(T) takes a value 
0.201(2). Remarkably, this value is significantly smaller 
than the value of the conventional Kosterlitz-Thouless 
(KT) transition 77 = 1/4 [7]. This observation suggests 
that the universality of the spin-transition of the present 
2D FFXY model might be different from that of the the 
conventional KT transition. 

To get further insight, we also examine the helicity 
modulus. In 2D XY models, the helicity modulus is 
expected to exhibit a discontinuous jump from zero at 
higher temperatures to a nonzero value T s / (2Tn](T s )) at 
the critical temperature T s . From the raw data dis- 
played in Fig. [T31 we extract the intersection points 
between the helicity-modulus data of each size and a 
straight line T / (2nr](T s )) by setting ?y(T s ) = 0.201(2) as 
obtained above. Then, the size dependence of the result- 
ing intersection points is displayed in Fig. 1141 Assum- 
ing that the correction term scales as l/ln 2 (L), we find 



8 



T 



0.465 
0.46 

0.455 
0.45 

0.445 
0.44 



- 
- 


§s( L 
|s(L 


v i 

4L) 
8L) 


1 ■— 



... .A.- - 




















■ 




- 


- 








..-*'" 


if'' 




• 










"t 











0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 

1/ln (L) 2 



FIG. 10: Crossing temperatures of the spin Binder pa- 
rameter g s (Fig. IHTa)). of the spin correlation-length ratio 
£s/L (Fig. [Hfb)), and of the effective vorticity modulus vi 
(Fig. 0(a)) are plotted versus (1/hiL) 2 . L in the horizontal 
axis is the average system size. The combined straight-line 
fits for three kinds of crossing temperatures assuming a com- 
mon T s yields T s = 0.4418(5). In addition, crossing temper- 
atures of the spin correlation- length ratio £ S /L between the 
two sizes L and sL with s = 2,4,8 are also shown as a func- 
tion of l/(lnL) 2 . The combined fit using the last three data 
sets assuming a common T s yields T s — 0.4414(3). 



T s = 0.44185(5). This value of T s is consistent with our 
estimate above T s = 0.4418(5) obtained from the Binder 
parameter, the correlation-length ratio and the vorticity 
modulus. 

If one assumes instead the conventional KT value 
r?(T s ) = 1/4 in the above analysis of the helicity mod- 
ulus as was done in Refs.fl6l. [2lL |22| . one obtains T s = 
0.44619(4) as displayed in Fig. Q3| The latter value 
agrees well with the estimate of T s = 0.446 (Ref.fl6|) 
and T s = 0.4461(1) (Ref. 0, S2] ) . However, this value is 
not compatible with our present estimate of T s obtained 
from other quantities without assuming a specific value 
for the exponent T](T B ), T s = 0.4415(3). 

Our result indicates that the jump of the helicity- 
modulus is greater than the so-called 'universal jump' 
for the conventional KT transition [39l - l4ll |. Lee and Lee 
[ll| made a similar observation in the FFXY model on 
the triangular lattice. Indeed, such a possibility was men- 
tioned occasionally since early days of research @ ■ 



IV. SUMMARY AND DISCUSSIONS 

We studied the ordering of the spin and the chirality 
of the 2D FFXY model on the square lattice by extensive 
Monte Carlo simulations performed up to very large sys- 
tem size to circumvent crossover effects. We found the 
the chiral transition takes place at T c = 0.45324(1), while 
spin transition takes place at T s = 0.4418(5). Thus, the 
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FIG. 11: The sublattice spin susceptibility divided by L 2 is 
plotted versus L on a log- log plot. The straight lines represent 
power laws L~ v with n = 1/4 (green broken line) and with 
n = 0.201 (red solid line). 
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FIG. 12: The temperature dependence of the effective 
spin exponent ?7eft(T) extracted from the fitting of the sub- 
lattice spin susceptibility Xs(T,L) to the form Xs(T,L) = 
6 ( T ) L 2-„ off (T) uging the data of L > 32 , > 64 and > 128. 



model certainly exhibits the spin-chirality decoupling, 
where T s lies below T c by about 2.6%. 

As pointed out by Hasenbusch et al [H|,[22|], in order to 
observe separation of the two transitions it is crucial to 
simulate large enough system sizes. In Fig. [15] we display 
the raw data of the correlation length of the chirality and 
of the spin. Evidently, the two correlation lengths are al- 
most equal at higher temperatures T > 0.468, behaving 
side by side. This is the spin-chirality coupling behav- 
ior expected in the high-temperature regime. At lower 
temperatures, the bulk chiral correlation length becomes 
definitely longer than the bulk spin correlation length. 
For instance, at T = 0.456 we find £ c ~ 65 while £ s ~ 45. 
In the temperature range T < 0.458, we get into the 
spin-chirality decoupling regime. At the chiral transi- 
tion temperature T c = 0.45324(1), the chiral correlation 
length diverges while the spin correlation length £, S (T C , L) 
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FIG. 13: The temperature dependence of the helicity modu- 
lus. The straight lines represent T / (2^77) with 77 = 0.201 and 
1/4. 
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FIG. 14: Intersection points of the helicity modulus ob- 
tained by assuming T](T S ) = 0.201 and 1/4 are plotted versus 
l/(hiL) 2 . The straight lines represent the linear fits of the 
data. 



stays at a finite value around 120 as shown in the inset 
of Fig. El 

The universality class of each transition was also stud- 
ied. The universality class of the chiral transition is well 
compatible with the Ising universality. On the other 
hand, the universality class of the spin transition is dif- 
ferent from the conventional KT transition. 

Possibility of the non-universal jump of the helicity 
modulus was pointed out on the general ground of a re- 
fined rcnormalization-group scheme by Minnhagen (42j 
and was argued to be relevant to the fully frustrated XY 
model (HI. However the true connection to the frustra- 
tion system was not necessarily clear. The spin and the 
chirality orderings of the fully frustrated XY model was 
also studied by Choi and Stroud by means of a renormal- 
ization group analysis of the coupled XY model which is 
believed to be essentially equivalent to the FFXY model 



FIG. 15: Direct comparison of the chiral and the spin corre- 
lation lengths for the sizes L = 128 and 256. The full/open 
symbols are the finite-size data of the chiral/spin correlation 
length, respectively. 



[33|. The possible non- universal jump of the helicity 
modulus was suggested [13, HH . The latter RG scheme 
[33l — f35[ | is quite appealing since it aimed at studying the 
ordering behaviours of both the chirality and the spin in 
a unified manner. However, the assumption underlying 
the derivation of the RG equations, i. e. the smallness of 
the fugacity of the both the usual spin vorticies and ex- 
tra charges responsible for the ordering of the chirality, 
each being justifiable in the low-temperature and in the 
high-temperature limits, cannot be satisfied asymptoti- 
cally upon rcnormalization at any temperature. It will be 
very interesting to investigate further the nature of the 
spin transition in the presence of the chiral long-ranged 
order, which is not considered in conventional theories 
for the KT transition. 

The numerical calculations were carried out partly on 
SX8 at Yukawa Institute of Theoretical Physics in Kyoto 
University and on SX-8R at the Cybermedia Center in 
Osaka University. This work is supported by Grant-Aid 
for Scientific Research on Priority Areas "Novel States of 
Matter Induced by Frustration" (19052006). 



Appendix A: Vorticity modulus in the coulomb gas 
representation 

In order to get some insight into the meaning of the 
vorticity modulus, we derive in this appendix its coulomb 
representation. 

The partition function of our model can be represented 
as, 



i <ij> 



-pE i3 ({6 i3 },{A XJ }) 



(Al) 



with 



E l3 ({0 ij },{A ij }) = -Jcos( 



Aij), (A2) 
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The parameter Ay (vector potential) is understood as an 
anti-symmetric tensor, 



A, 



-A 



(A3) 



For a fcrromagnctic/antifcrromagnctic bond, we choose 
Aij =0/tt, respectively. 

We consider the change in the free energy by an in- 
finitesimal increase of the vorticity m, 



V 



d 2 F{m) 



dm 2 



(A4) 



m=0 



where F(m) is the free-energy of the system with a mod- 
ified vector potential, 



A, 



A> 



(A5) 



where <pij is an 'angle' between the sites i and j with 
respect to the origin O, as shown in Fig. [16] We set, 

<t>ij = ~<t>ji- ( A6 ) 




FIG. 16: Solid angle </>y between sites i and j with respect to 
the origin O. 

By the standard steps of mappings [40j ] . we find that 
the partition function of the system can be represented 
as, 



Z = Z sw ]~ e" 



(A7) 



where Zgw is the contribution from spin-waves and 
H ({qn}) is an effective charge Hamiltonian for the charge 
variables q^ defined on the dual lattice with its lattice 
points n = (x n ,y n ) located at the centers of the plaquc- 
ttes of the original lattice, 

qn = qa + bn- (A8) 
Here q^ is an integer and bfi is a 'flux' defined by 



— y a 

<ij> around n 



(A9) 



where the sum is taken around the plaqucttc at n in a 
counter-clockwise direction. 



One can see that the contribution from the additional 
term m4>ij in the vector potential Eq. (|A5[) to the flux ha 
is zero almost everywhere except at the origin ft = O; 

bH = m6 H o + — ]T A i:j . (A10) 

<ij> around n 

Here it is made evident that m can be interpreted indeed 
as an infinitesimal charge or vorticity. 
We then find, 



d 2 (-f3F) 



dm 2 



m =0 \ 9(q 5 ) 2 / 



where 



d{-m) \ 



-mm) 



1 In t—^ Qn = — c 



-mm) 



(All) 



(A12) 



An explicit form of the charge Hamiltonian 1-L{{q}) can 
be written as, 



^({9})= 92{ni,n 2 )q i t 1 qn 2 



(A13) 



where higher-body (i.e., 4-, 6-,. .. body) terms are ne- 
glected for simplicity, which corresponds to the Villain's 
approximation. The two-body interaction 32(^1,^2) is 
essentially the same as the spin-wave Green's function, 

<? 2 (ni,n 2 ) ee E 4 _2( C os(fc x )+cos(fc y )) 



k 
1 

2^ 



ik\m — n 2 | 



dkk- 



/L 



(A14) 



In the first equation, the summation over all possible 
wavevectors k — (k x ,k y ) allowed in the periodic system 
of size L x L is taken. The Green's function behaves at 
large distances as, 

5 2 (ni,n 2 ) «27rln(|n 1 -n 2 |)+2 7 rln(2V2e 7 ), (A15) 

where 7 is Eulcr's constant while 

g 2 (n,n) w 2tt ln(i) + 2tt ln(2 v / 2e 7 ), (A16) 

which is the bare vortex core energy of a finite sized sys- 
tem of size L. 

Within the Villain's approximation, we find, 



V 



d 2 F 



dm 2 



32(0, 5) - P n)g 2 (0, n')({q ri q-,)c - (fe)c(^)c)- 



(A17) 
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The first term on the r.h.s of the last equation is obvi- 
ously the vortex creation energy which generally scales 
as c(0) + v(0) ln(i) with a 'bare' vorticity modulus at 
zero temperature v(0) = 2ir and a constant c(0) = 
2tt ln(2v / 2e 7 ). The second term would be interpreted as 
a 'screening' term due to thermal fluctuations of other 
charges. On general grounds, we speculate that the same 
scaling holds at finite temperatures below the critical 
temperature T s , 

V = c(T)+v(T)]n(L), T < T s , (A18) 



with a temperature-dependent vorticity modulus v(T). 
Above T s , we have v{T) = in the thermodynamic limit 
L — > oo. 

The obtained expression of V, Eq. (| A17I) . should be 
compared with the correspondin g ex pression for the he- 
licity modulus which is given by [161 |4l| , 

T = 1 - ^ 9%% n'Wnq^c - fe)cfe)c). (A19) 

ft,n' 
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